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THEORETICAL AERODYNAMICS OF UPPER-SURFACE -BLOWING 

JET -WING INTERACTION 

C. Edward Lan* and James F. Campbell 
Langley Research Center 

SUMMARY 

A linear, inviscid subsonic compressible flow theory is formulated to treat the 
aerodynamic interaction between the wing and an inviscid upper-surface-blowing (USB) 
thick jet with Mach number nonuniformity. The predicted results showed reasonably good 
agreement with some available lift and induced-drag data. It was also shown that the 
thin- jet -flap theory is inadequate for the USB configurations with thick jet. 

Additional theoretical results showed that the lift and induced drag were reduced 
by increasing jet temperature and increased by increasing jet Mach number. Reducing 
jet aspect ratio, while holding jet area constant, caused reductions in lift, induced drag, 
and pitching moment at a given ai^le of attack but with a minimal charge in the curve of 
lift coefficient against induced-drag coefficient. The jet -deflection effect was shown to be 
beneficial to cruise performance. The aerodynamic center was shifted forward by addii^ 
power or jet -deflection angle. Moving the jet away from the wing surface resulted in 
rapid chaises in lift and induced drag. Reducing the wing span of a rectangular wing by 
half decreased the jet-circulation lift by only 24 percent at a thrust coefficient of 2. 

INTRODUCTION 

Recently, the feasibility of using the upper-surface-blowing (USB) concept for short 
take-off and landing operations has been under extensive experimental investigation (refs. 1 
to 7). Some of these experimental results have been summarized in a recent paper by 
Johnson and Phelps (ref. 8), These experimental investigations confirmed the promisir^ 
aerodynamic and noise characteristics of the USB configurations at low speeds under high 
lift conditions. On the other hand, several studies, such as references 9 and 10, have 
indicated that design for good low-speed performance, for example, spanwise spreadii^ 
into a thinner jet for good jet turning, would appreciably compromise the cruise perfor- 
mance. Additional information about the jet influence is needed, therefore, to provide the 
designer more choices in his trade-off study. Although the USB concept has been in exis- 
tence for quite a while , only a few limited theoretical methods have been developed to 
describe the associated aerodynamic phenomena. 

^Associate Professor in the Aerospace Ei^ineering Department at the University of 
Kansas, Lawrence, Kansas. 



Experimental evidence shows that the USB jet flow is much more concentrated 
than that associated with the external-blown-flap (EBF) configuration. It is well known 
that the conventional thin jet flap is much more effective in producii^ the lift for a given 
total jet momentum if the blowing jet has been spread out over a large spanwise distance. 
This fact may explain some success in predicting the lift characteristics of the EBF con- 
figuration by thin-jet-flap theory (ref, 11). The high lift capability of the USB configura- 
tion, however, can hardly be explained with the thin-jet-flap theory, because the jet span 
of the USB configuration tends to be relatively narrow and the jet is thick. Some calcu- 
lations using the thin-jet-flap theory with the narrow jet span showed that it consistently 
imderestimated the measured lift for the USB configurations. 

Original USB concepts, similar to the one in reference 12, had jet thicknesses at 
the jet exit that were relatively small, less than 2 percent of the local chord. In recent 
experiments jet thicknesses were examined which were greater than 10 percent of the 
local chord, a result of simulating the modern or future high-bypass -ratio turbofan 
ermines. With a jet-flow region of finite thickness and of higher dynamic pressure than 
the free stream, it is not unusual that additional lift would be induced on the neighboring 
lifting surfaces. This fact has long been recognized by Kuchemann and Weber (ref. 13). 
Physically, the wing-induced flow field will be modified by the jet region and the jet flow 
will be perturbed because of the presence of the wing. It is the neglect of this wing-jet 
interaction process which may account for the imderprediction of lift by the thin-jet-flap 
theory. 

Historically, the interaction between two flow fields with different energy levels 
has been concerned with the wing- propeller slipstream interaction problems. Some suc- 
cess in numerical computation has been reported. (See, for example, refs. 14 and 15.) 
Except for the study of reference 15, the Mach number effect has usually been neglected. 
Since in the wing-jet interaction problem, as being considered here, the jet Mach number 
can be different from the free-stream value (Mach number nonuniformity), the theoretical 
formulations similar to that of reference 14 cannot be conveniently used. Instead, the 
formulation of reference 15 will be adopted here. 

The purposes of this paper are, therefore: (1) to present a thick- jet- wing inter- 
action theory which accounts for differences between the jet and free-stream Mach num- 
bers, (2) to apply this theory to predicting the aerodynamic characteristics of USB config- 
urations, and (3) to investigate theoretically the effects of some configuration and jet 
parameters. Linear, inviscid subsonic compressible flow theory will be assumed. The 
effects of wing thickness, nacelle, and fuselage have not been included in the analysis. 
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SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements and cal- 
culations were made in U.S. Customary Units. 


m 

wing aspect ratio 



jet aspect ratio 



jet cross-sectional area, m^ (ft^) 


a 

= (x^ - x)T + (y^ - y)T + (zi - z)k, m 

(ft) 

a' 

= (x^ - x)i + ^(yj - y)j + - z)k, 

m (ft) 

b 

span, m (ft) 


b 

= (xg - x)i + (y2 - y)T + (zg - z)k, m 

(ft) 

b’ 

= (x2 - x)T + ^(y2 - yjT + ^(Z2 - z)k. 

m (ft) 

C 

leading-edge suction parameter (see 

eq. (41)) 

^D.i 

induced-drag coefficient 


^D,j 

drag coefficient due to jet deflection 



Cl total lift coefficient (circulation lift plus jet-reaction lift) 

-^^Cl difference in lift coefficients with jet on and off, except as noted 


^L,r 

^l,T 

Cm 


circulation lift coefficient 
circulation rolling-moment coefficient 
pitching-moment coefficient about Y-axis 


Cm,r pitchii^ -moment coefficient about Y-axis due to circulation only 
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M,i 

Cf 

‘^jet 


'£,r 


-m 


’m,r 


D 


L 


net thrust coefficient of jet, Net thrust/q^S 
jet-momentum coefficient, mVj^q^S 
local chord length, m (ft) 
reference chord length, m (ft) 
sectional induced-drag coefficient 
flap -chord lei^th, m (ft) 

local chord length at jet center line (see fig. 1), m (ft) 

sectional lift coefficient 

sectional circulation lift coefficient 

sectional pitching -moment coefficient about Y-axis 

sectional pitching-moment coefficient about Y-axis due to circulation 
only 

sectional leading-edge thrust coefficient 
sectional jet-momentum coefficient 
drag, N (lb) 

unit vector tangent to jet path (see fig, 1) 
vector defined by equation (A1 7) 

distance of jet lower surface to wing surface, m (ft) 
unit vectors along X-, Y-, and Z-axes, respectively 
total lift, N (lb) 
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= “ *1? + (^2 " yi)5' + i^2 - 

= (xg - xi)T + ^(y 2 - yi)T + ^(zg - Zj)k, m (ft) 
Mach nvimber, or number of integration points 
vector defined by equation (A3) of appendix 
mass flow rate, N/sec (slv^s/sec) 
nxunber of chordwise integration points 
normal velocity influence -coefficient matrix 
unit vector normal to jet surface (see fig. 1) 
jet axis system (see fig. 1) 
static pressure, N/m^ (Ib/ft^) 
dynamic pressure, N/m^ (Ib/ft^) 
induced-velocity vector, m/sec (ft/sec) 
radius of curvature, m (ft) 

= xi + yj + zk, m (ft) 

= (x - x*)^ + ^2(y - y’)^ + )S2(z - z*)^, m^ (ft2) 

o n 

wir^ area, m (ft ) 

tangential velocity influence-coefficient matrix 
jet -path coordinate, m (ft) 

= Po/pj 

jet thickness, m (ft) 
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nondimensional perturbed velocity in X-direction 


velocity 

unperturbed velocity vector, m/ sec (ft/ sec) 

perturbed velocity vector, m/sec (ft/sec) 

jet-entrained -flow velocity vector, m/sec (ft/sec) 

wing -fixed rectangular coordinates with positive X-axis along axis of 
symmetry pointing downstream, positive Y-axis pointing to right, 
and positive Z-axis pointing upward, m (ft) 

coordinate of camber surface, Zg = Zg(x,y), m (ft) 

angle of attack, deg 


nondimensional vortex density, or ratio of specific heats 

flap deflection, deg 

jet -deflection angle, deg 

angular coordinate (see eq. (A5) of appendix) 

angle of camber slope 

sweep angle, deg 

reflection coefficient as defined in reference 24 



p 

4 > 

T 

4 ' 

Subscripts: 

a 

e 

i 

jj 

LE 

0 
oj 
TE 
w 
wa 
wj 
wo 

r 

1 


density, kg/m^ (sl\igs/ft^) 

dimensional perturbation velocity potential, m^/ sec (ft^/ sec) 
nondimensional perturbation velocity potential 
nondimensional additional perturbation velocity potential 

additional 
jet exit 
jet flow 

jet vortices for jet flow 
leading edge 
outer flow 

jet vortices for outer flow 

trailing edge 

wiig 

additional wii^ vortices 
perturbation due to wing in jet flow 
perturbation due to wing in outer flow 
circulation 

first end point of vortex element 
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second end point of vortex element 


oo free stream 

MATHEMATICAL FORMULATION 


Boundary Conditions 

The unperturbed flow, that is, the flow field before the wing is introduced, includes 
in general the uniform free stream, the jet flow, and the jet -entrained flow due to turbu- 
lent mixing. In the following formulation the jet-entrained flow may be regarded as 
being combined with the uniform free stream to form the outer flow represented by the 
vector Vq. As the wing is introduced into the unperturbed flow, perturbations of the jet- 
flow and the outer-flow fields occur. The calculation of these perturbations represents 
one of the main problems in the present formulation. The unperturbed jet flow is 
assumed to be uniform. 


These perturbations can be calculated by satisfying the boundary conditions on 
both the jet and the wing surfaces. The jet boundary conditions require that the jet sur- 
face be a stream surface and the pressure be continuous, that is, the pressures on either 
side of the surface be the same. The stream surface condition is satisfied if the slopes 
of streamlines on both sides of the jet surface are the same, that is. 


n • (Vq + Vq) _ n • (Vj + Vj) 
e-(Vo+ Vq) e .(Vj +fj) 


where e is a xmit vector tangent to the jet path and n is a unit normal vector at the jet 
surface. (See fig, 1.) The unperturbed jet flow is represented by Vj and the perturba- 
tions are represented by small letters with appropriate subscripts. 


The pressure continuity condition can be formulated (see ref. 15) by using 
Bernoulli's equation so that 

y. 


1 


2 ’'jpj 


V- + V- 




= Pr 


1 + 


y - 1 Pq 
2 yp^ 


2 


y 

-ly-i 


Vo + Vo 


( 2 ) 


It is assumed that in the unperturbed flow the jet static pressure, Pj, at the jet exit is 
equal to the pressure of the outer flow, p^, at the jet exit, so that the static pressures 
are matched. In addition, if yj = y, equation (2) can be reduced to (see ref. 15) 
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(3) 



+ 





V. = P + -5- 


To simplify equations (1) and (3), it may be assumed that the flow fields are irro- 
tational. Let 

Vj = (Yq - 11 ) 11 + V<^j (4) 

Vo = (5) 

where ^Vq - n) is the component of the unperturbed outer flow normal to the jet surface. 
Introduction of equations (4) and (5) into equation (1) yields 




assuming that e • V0j « e • Vj and • V<^^« e • where 


.•-JjLll'' 


V 


i'® 


Vj-n = 0 


( 6 ) 


(7) 


Equation (6) is the linearized flow tai^ency condition at the jet surface. It is seen to be 
the same as equation (8) of reference 16 in wing- slipstream interaction formulation. 

Similarly, a linearized version of the pressure continuity condition is (see ref. 15) 


P 



as 




( 8 ) 


where s is the distance measured along the jet path and is taken to be x in the shallow- 
jet assumption. In derivir^ equation (8), the quantities . n and Vj • n are assximed 
small. 

The wing- flow tangency condition will now be considered. The wing is assumed 
to be completely in the outer-flow region only. The conventional flow tangency condition 
is therefore applicable. It follows that on the wing surface the following condition must be 
satisfied: 


dch 

dZ 



where z^(x,y) is the wing-camber fimction. 


(9) 
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Method of Flow Singularities 


To satisfy the boundary conditions, equations (6), (8), and (9), the method of flow 
singularities will be used. For this purpose, it is convenient to nondimensionalize 
the boundary conditions. Let the nondimensional perturbation potentials (p^ and 
be defined such that 



(10) 


<p.= Y..e U> 


( 11 ) 


If it is assumed that 


9(Vo . e 


and 


• e) 


9n 9s 

tions, equations (6), (8), and (9), become, respectively. 


are small or zero, the boundary condi- 


where 


9n 


9n 


Vo • n(l - /,’) 



9s 


t(m’) 


2 

9s 


= 0 


_ Vq • 1 9Zg ^ Vq ’ h 
9z Vq . e 9x Vq . e 



(On jet surface) 

(12) 

(On jet surface) 

(13) 

(On wing surface) 

(14) 


(15) 


For numerical convenience, it is advisable to evaluate the wing-alone case sep- 
arately. For this purpose, it may be noted that in equation (14) Vq • T « Vq ' e in the 
present linear formulation. Let Vjg be the jet-entrained flow and V^ the uniform 
free-stream vector. Then 


Vq = V^ + Vje (16) 

It is assumed that Vq • e" « V^ cos a. 

Two nondimensional potential functions are defined for the wing -alone situation, 
where there is no jet flow, and are illustrated in the following sketch: 
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Region occupied by jet 



<P. 


wj 



0 

wo 


y 

4^ 


/ 


Wing 


The symbol represents the wing in a uniform free-stream flow having Mach 

number Mq. The symbol is for the wing in a free-stream flow having Mach 

number Mj and is defined in that region of the wing's flow field to be occupied by the 
jet. The presence of the jet is accoimted for by the additional perturbation potentials, 
and , which are specified for the regions outside and inside the jet. (See the fol- 
lowing sketch.) 



M, 






L 


Wing 


These additional potentials are added to the potentials for the wing-alone flow fields to 
obtain the total potentials inside and outside of the jet region. Thus, 


^ - <’wo(“o) + (11) 


Using equations (17) and (18) the boundary conditions, equations (12) to (14), become 


an 


an 


Vq . n(l - fi’) 


Vo*^ 


an 



Q<^wo 

an 



(19) 


11 



( 20 ) 


8^. 

as 


T(fx’)^ 


9s "* 9s 


T(ju’)^ 


d(b 

^wo 

“a^ 



az cos a 

OO 


(On wing surface) 


( 21 ) 


where and both satisfy the condition that the flow be tai^ent to the wii^ cam- 

ber surface: 


^‘^wo(^o) 

az 


azc ^ 

- tan a 
ax 



9zc 

ax 


- tan a 


(22a) 

(22b) 


The formulation given in equations (19) to (21) is convenient in that, if ju' = 1, = Mj, 

and V. = 0, then the additional perturbation potentials will be automatically zero. 

At this point, a two- vortex- sheet model for the jet surface will be introduced, as 
illustrated in the following sketch: 







r> o 




\jf ^ 

J2L - 


One sheet is to account for the additional perturbation in the outer -flow field and the other 
for the jet -flow perturbation. This arrangement is required in order to represent the 
flow field correctly. Any singularity distribution should induce different normal and tan- 
gential velocities at all other points on the jet surface depending on whether the surface 
is approached from inside or outside the jet region, due to Mach number nonuniformity. 
Ribner and Ellis (ref. 14) introduced a source distribution to account for the jump in nor- 
mal velocity (eq. (19)). This source distribution would be dependent on the vortex distri- 
bution introduced to satisfy equation (20) and vice versa. Hence, additional new integrals 
for the source distribution must be handled. On the other hand, the present method calls 
for the use of two vortex sheets so that only vortex integrals must be dealt with. 

The resulting vortex integrals are now reduced to finite sums through a quasi- 
vortex-lattice method (ref. 17). This method represents a refined vortex- lattice method 
by accounting for the wing-edge singularities and Cauchy singxilarity in the integrand. 

It has been found to be more accurate than the conventional vortex- lattice method for 
a given number of vortex elements used. This property of faster rate of convergence 
of the method is important in the present formulation in that without it the number of 
unknown vortex elements would be greatly increased. Now [Sj^ may be defined as 
the s- influence-coefficient matrix for the jet tangential velocities due to the wing; that is. 
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the matrix elements being the s- perturbed velocities on the jet surface due to unit vortex 
elements on the wing. Similar definitions are applicable to the matrices 
[Nw^, etc., where N denotes the normal components. Using y with appropriate 
subscripts to denote the unknown vortex strength, additional perturbed velocities on the 
jet surface can be written as 



I^Jw](o)(^wa) + |?Jj]^^^{’^oj} 

(23a) 


(23b) 

[®Jw](o)0^wa) + &jj]^^^{^oj) 

(23c) 

' C/3 ' 
! 

(23d) 


where the subscripts o and j denote the outer and the jet regions, respectively, 
and represents the additional wing vortex strengths induced by the jet. Simi- 

larly, on the wing surface. 


= [Nww](>^wa) + [^Wj], .{^oj} 


(o)' 

Substitution of equations (23) into equations (19) to (21) gives 

Vo • n(l - (i') . 9^wj 


(23e) 




30 


(o) 


*(j)^ 


Vo-e 


9n 


(“)) - ^(“o) 




wo\ 


'(o) 


'(j)' 


3s 


and 


[^Ww]{>^wa) [%j]^^^{'^oj) - f 


cos a I 


(24) 

(25) 

(26) 


Equations (24) to (26) are written in the form of an augmented matrix equation 
and are solved by using Purcell’s vector method for solving simultaneous equations 


(ref. 18). The main advantage of this solution technique is that it requires only 
memory locations for the matrix operation, where is the matrix size. 


(Nmj^ 
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Once the additional wing vortex strengths, have been determined, the total 

wing y^ is, then. 


y = 

'Vf 


y + 
wo 


wa 


( 27 ) 


where y^^ is the wing vortex strength in the uniform flow. 


Effect of Jet -Sheet Curvature 

That part of the jet which is directly on the wing surface is assumed to have no 
curvature and is therefore independent of the flap deflection. It follows that, if the jet- 
entrained flow is neglected, the flow tangency condition for the jet becomes 


Vq ♦ n(l - 

Vo-e 


-(1 - fx’) tan a 


(28) 


The deflected jet sheet behind the trailir^ edge will be subject to varyii^ angle of attack 
due to the jet-sheet curvature. This curvature distribution is unknown and is dependent 
on the unknown jet vortex strei^h. To avoid iteration, Spence’s method in the thin-jet- 
flap theory (ref. 19) can be adopted here in the followii^ way. Since the jet flow is 
assumed irrotational, it is true that (see sketch) 



Vj + Vo 

where V. = -r . The subscripts 1 and 2 used in this section refer to conditions 

J ^ 

in the jet region on the upper and lower jet surfaces, respectively. Introducing the non- 
dimensional perturbed velocity potential as given in equation (11), it can be seen that 


V 


1 


= V. + V. 
1 ] 



V 


2 


^3 ^ 
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and hence, 



Substituting equation (13) into equation (30) yields 


R 




= T(m’)^ 


. 9s/ 


^O 

V 9S 


f B(j) 


WO 


BS 


80. 


WO 


BS 


For convenience, let the right-hand side of equation (31) be denoted by f(x): 


( 30 ) 


(31) 


f(x) = T(/x’)^ 





(32) 


1 /i2 

For small jet deflection, « 5-|. It follows that the jet path satisfies the followii^ 

K (jx^ 

initial- value problem: 





(33) 


where 5j is the jet -deflection angle. Since only the slope ^ is needed, complete inte- 
gration of equations (33) for z(x) is not necessary. To integrate, let 


X 


Cl 

Xte + y(l - cos e) 


(34) 


where Cj is the jet length from the wii^ trailing edge to be included in the analysis. 
Then equations (33) become 
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Hence, 




^0 

A9c4 






/ y sin 0^1! 

x(0k) 

+ i sin 0^f 

- 1 




1 


where the integral has been reduced to a finite sum through the trapezoidal rule. It fol- 
lows that at the ith control point on the trailing jet sheet, the jet-path slope is given by 


dx/j^ 


A0C,- 




2t 


^ k= 1 


sin 0}jf 


X(0T 


+ ~ sin 


:^9l) 


(35) 


Once ^ for the jet path is obtained in terms of the unknown vortex strength, 
the first term on the right-hand side of equation (24) can be written as (again neglecting 
the jet-entrained flow) 


V„ • n(l - p’) 


dz 

dx 


tan cvl'^1 - 




(36) 


j _ 

Those terms in ^ which depend on the unknown vortex strength are then combined into 
the left-hand side of equation (24) before the solution is attempted. 


Since the s- induced velocities are needed in equation (32), and therefore in the 
jet- surface tangency condition, equation (24), it is convenient to solve equation (25) first 
and then equations (24) and (26). This procedure can also avoid the numerical difficulty 
discussed in reference 15 when some of the jet-control points coincide with some of the 
wing-control points. The latter would happen if the jet is blowing on the wing surface. 


Calculation of Aerodynamic Forces and Moments 

To find the sectional aerodynamic characteristics, it is assumed that the wing 
vortex elements are situated along the actual camber surface. Since the resulting pres- 
sure force is normal to the camber surface, the sectional aerodynamic characteristics 
can be determined by resolving the pressure force in the proper direction and integrating 
over the local chord. It follows that 
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where 


""a = r ^oVo(^o • + ‘'t (“ - ^le) 

p V c 
'^o 


^ 2cos«y^ 

- 0(x) dx + sin^O! - ^le) 

(37) 

%i= Vx)sin[ 

01 - 0(x^ dx - cos - ^le) 

(38) 

c * 2_gos a r y (x)x cos 
cc ^0 ^ 

[a - ^(x^ dx 

(39) 


9(x) = tan 



(40) 


The integrals in equations (37) to (39) are reduced to finite sums through a modified trap- 
ezoidal rule after they are reduced to an integration over an angular coordinate (see 
ref. 17). 

The sectional leading-edge thrust coefficient Cj. is computed as (ref. 17) 




1 1 \/i 

2 cos Al 2 V 



cos 



(41) 


where C is the leading-edge suction parameter. According to reference 17 the param- 
eter C can be determined by smnming the total induced normal velocity at the leading 
edge and subtracting the right-hand side of equation (22a) or equation (26), as the case 
may be. If the induced downwash at the ith leading- edge control point is denoted by a sub- 
script i, then in a manner similar to that used in reference 17, the uniform flow case is 
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Similarly, for the additional flow, it can be shown from equation (26) that 


NC2 



^LE 





{ 


V cos a I 

oo ' 


( 43 ) 


The suction parameter is then C = + C 2 . 

The overall aerodynamic characteristics of the wing are determined by spanwise 
integration of the sectional characteristics, as described in reference 17. 

The expressions for the elements of the influence-coefficient matrices are described 
in the appendix. 


Numerical Considerations 

Since the u-induced velocities, or the [s] matrices, are needed in the present for- 
mulation, special care must be exercised in their evaluation. To illustrate this point, 
consider the expression for u(x,z) for the two-dimensional case in incompressible flow: 



y(x') dx' 
(x - x')2 + 


(44) 


Observe that for 0 <x < 1 the integrand in equation (44) has a second-order singularity 
as z approaches zero while the whole integral is multiplied by z. In fact, it can be 

shown that u(x,0) = As has been checked numerically, equation (44) can be inte- 

grated accurately by the integration technique in the quasi-vortex-lattice method (ref. 17) 
only if z and the number of integration points are not too small. Increasing the number 
of integration points, however, would become numerically unrealistic because the number 
of integration points is equal to the number of unknowns to be solved. 

As has been shown in reference 17, the integration technique of the quasi-vortex- 
lattice method is quite efficient in treating integrals with Cauchy singularity. The best 
numerical technique under the present circumstances, therefore, is to rewrite equa- 
tion (44) as 


u(x,z) = I- r' dx' + f ' 


dx’ 


^ ^ 2L Y 
^ 47t N /j 


(x - x')*^ + z^ 
N 


rfSk) " 


2^^ ^0 (x - x’)2 + z2 

M 


sin 0. + 


k= 1 (x - Xj^) + z^ 


JL y 

M Zy 


sin 0j 


zy(x) 

4" MA,(x-Xjf.z2 


(45) 
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where 


9, = 2k - 1 ^ 
k 2N ^ 


Xk = |(l - cos 0k' 




XJ = i(l - cos 9j) 


and M must be chosen so that 


(k » 1, 2, . . N) 

i (46) 

(j « 1, 2, . . M) 


M = 2^N (47) 

for interdigitation between the control and integration points in the last summation and p 
is an arbitrary integer. This procedure has been found to be quite accurate in most prac- 
tical applications (ref, 15). 


RESULTS AND DISCUSSION 

In order to compare the present aerodynamic theory with some available experimen- 
tal data, it is necessary to account for the reaction forces associated with the jet. The 
following sketch illustrates the ram drag, at the inlet of the nacelle and the jet- 

reaction force, niVj, at the trailing edge of the wing. 



When the flap is deflected, the jet will also deflect due to the viscous interaction 
between the jet and the flap. Since this effect, commonly referred to as the Coanda effect, 
can be explained only by viscous -flow theory, its exact prediction is beyond the scope of 
this study. For the present calculations, therefore, the jet is assumed to deflect with the 
flap and to preserve its momentum to the trailir^ edge. This assumption results in a jet- 
reaction force which is combined with the ram drag to obtain the following jet-deflection 
forces in the lift and drag directions: 

Lj = mVj sin ^6^ + aj (48) 
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(49) 


+ aj 

In the following applications, the jet is assumed to have a constant width, thick- 
ness, and velocity, and the jet -entrained flow has been ignored. The lateral expansion of 
the jet with wind on has been observed experimentally to be relatively limited (refs. 8 
and 9). Even if the jet lateral expansion is not ignored in the theoretical model, the 
increase in lift due to the larger washed area of the lifting surface will be compensated 
by the associated decrease in the jet dynamic pressure due to enlargement of the jet 
cross section. A similar situation has been observed experimentally by Seidel (ref. 20) 
using circular jets with the effect that the wing lift is nearly independent of the horizontal 
nozzle distance. Similarly, if the jet is allowed to expand vertically, the increase in lift 
due to larger jet thickness will also be compensated by the associated decrease in the jet 
dynamic pressure. (See refs. 4 and 21.) It should be emphasized that these are simply 
assumptions in the present applications of the theory and have not been vigorously verified 


Dj = riiV^ - mVj cos ^5j 


Comparison With Experimental and Theoretical Thin- Jet- Flap Results 

There are only limited experimental data for USB configurations available for veri- 
fication of the present theoretical model. Extensive theoretical and experimental results 
are available, however, for jet -flap configurations, where it is assumed that a thin jet 
sheet emanates from the trailing edge. The computer program was slightly modified, 
therefore, to treat the thin-jet-flap case. According to the thin-jet-flap theory, it is 

Vi 9 

assumed that — L oo but V. t. is finite, where t^ is the jet thickness. It follows that 

Vqq j j 

jU = 0 and the jet flow is imperturbable. This fact implies that — ^ = 0 in equation (12). 

9n 

In equations (31), (32), and (35), furthermore, the expression can be simplified 

as 


T(|i')^ 2 

p.v?t. c^(y)c 


(50) 


PjVfti 

where c^(y) = ^ ^ the sectional jet -momentum coefficient. To avoid difficulty 

2 ^ 00^00 ^ 

with the presence of the jet sidewall vortices, the jet thickness was assumed to be 5 per- 
cent of the chord length. The trailing jet length to be included in the analysis was deter- 
mined numerically by using various lei^ths of the jet sheet until a maximum value for the 
predicted Cl was reached for a given number of streamwise vortices. From this 
nvimerical experimentation, it was found that for practical purposes a jet distance of 2 to 
3 chord ler^ths and 5 to 7 streamwise vortices are sufficient for numerical convergence. 
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Das' experimental and theoretical section lift and pitching- moment distributions for 
a wing with half semispan blowing and 6j = 30° (ref. 22) are compared with the present 
method in figure 2. Experimental and theoretical lift increments due to jet blowing for a 
rectangular wing with a partial-span flap (ref. 23) are compared with the present method 
in figure 3, The results by the present theory presented in figure 3 were obtained by 
adding the predicted incremental lift to the lift at = 0 given in reference 23. The 
comparisons in figures 2 and 3 indicate that the predictions of c^ and Cj^^ made with 
the present method are in good agreement with existing thin-jet-flap theories and exper- 
imental data. 


Comparison With Thick-Jet Experimental Results 

The interaction between a propeller slipstream and a wing is one type of thick-jet 
problem that has been examined using the present theoretical formulation. This compar- 
ison was reported in reference 15, which demonstrated that the current theory provides 
good agreement with existing theories and experimental data. In order to evaluate the 
predictive capabilities of the thick-jet-wing interaction theory for USB configurations, 
the present method is compared with the experimental results of references 4 and 1 in 
figures 4 and 5, respectively. In order to make these comparisons, the following condi- 
tions were imposed: 

(a) The jet -deflection angle was assumed to be equal to the flap angle, xmless other- 
wise noted. 


(b) The net thrust was assumed to be equal to the experimentally determined static 
gross thrust. 

(c) The jet-exhaust velocity was computed by using the momentum principle which 
neglects any friction losses, wake rotation, etc. In addition, the jet-exit static pressure 
was assumed to always be equal to the free-stream static pressure. The momentum 
theory leads to the followii^ expression which was used to compute velocity as a function 
of the net thrust coefficient: 



1 + 


2Ct(S/2)^ 

^i(Pj/Po) 


nl/2 


(51) 


where Aj is the jet cross-sectional area. (^Note that Crp is the total thrust for two 
engines. ) 


According to the experimental observations of references 4 and 9, the jet is deflected 
even though the flap is undeflected, mainly due to the curvature of the upper surface of the 
airfoil. The trailing- edge angle of the airfoil's upper surface was quoted in reference 4 
to be 12°; this was also the jet-deflection angle assumed in the present computation. The 
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configuration which has twist and camber is idealized as shown in figure 4(a). The theo- 
retical results at two Cx's compared with the measurements and thin- jet- flap results 
are illustrated in figure ’4(b). The experimental results are given in figure 36 of ref- 
erence 4. The theoretical curves are obtained by adding the predicted lift increments 
at each a to the experimental wing-alone results. The present method slightly overpre 
dieted the lift at Crj, = 2,0 and underpredicted it at C-p = 0.5. It is seen that the thin- 
jet-flap theory imderpredicted the lift for C-j- = 2.0. 

Since in the experimental data the thrust and scrubbing drag have all been included 
in the total drag, it is not possible to compare the induced-drag prediction alone. The 
following method can be used to check the theory indirectly, however. If it is assumed 
that the scrubbing drag is independent of a, then the following drag increment. 


AC 


D 


(tv) 


= C 


D 


- C 


{a) °( q '= 1 °) 


(52) 


can be regarded as due to the induced drag and the jet-reaction effect. The induced drag 
is directly computed in the program and the jet-reaction effect is estimated by equa- 
tion (49). The drag increment due to angle of attack, therefore, is estimated by the pres 
ent theory using the following relation: 


' (^D,i * - (CD,i + 

Using this procedure, the drag increment is obtained at a = 6° and 11° and is com- 
pared in table 1. It is seen that the agreement is reasonable. 

TABLE 1.- COMPARISON OF PREDICTED INDUCED DRAG INCREMENTS AT 

ANGLES OF ATTACK WITH MEASUREMENTS FROM REFERENCE 4 




^CD(q,) 


deg 

Cx == 0*5 

Cx = 2,0 


Theory 

Experiment 

Theory 

Experiment 

6 

0.043 

«0 

0.112 

0.15 

11 

.114 

.1 

.29 

.3 


The second test case is a configuration with full-span flap deflection of 30°, as 
given in reference 1. The idealized planform is illustrated in figure 5(a), where the body 
cutout, the leading- edge Krueger flap, and the trailing- edge flap extension have been 
ignored, and the jet-exit location was at 0.23 of the local chord. The lift-coefficient 
against angle-of-attack curves are compared in figure 5(b), where the theoretical curve 
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is obtained by addit^ the predicted lift increments to the experimental curve for power 
off. As with the first test case, the prediction of Cl is reasonably good. 

Configuration and Jet-Parameter Effects 

Effects of jet-exit l ongitudinal location. - The effects of moving the jet exit in the 
chordwise direction are illustrated in figure 6 for the configuration of figure 5(a). The 
figure shows the jet-induced lift generated by the jet, which is located at chord stations 
from the leading edge to the trailing edge, for three angles of attack. The results indi- 
cate that moving the jet exit forward from the trailing edge is beneficial for lift augmen- 
tation and that the largest lift augmentations occur with the jet exit located at, or forward 
of, the leading edge. 

Jet -thickness effects at constant C-p and jet span.- For all of the following calcu- 
lations, the transport- type planform shown in figure 7 will be used, unless otherwise noted. 
Assiuning a flat-plate airfoil and C^ = 0.2, to simulate a cruise condition, three sets of 
theoretical curves have been plotted in figure 8 to demonstrate the effects of jet thickness. 
The first and second curves allow the jet Mach number to vary according to 



It follows from equation (51) that, as the jet thickness is reduced with C^p constant and 
the jet width is fixed, the jet velocity and the jet Mach number would be increased. The 
third curve was obtained by takitg M^ = Mj = 0.3 and = p^. This assumption 
neglects the interaction process that takes place when the Mach numbers are different. 

The following observations can be made about the results shown in figure 8: 

(a) With the Mach number interaction between the streams neglected ^M^ = Mq^, the 
third curve shows that a jet thickness exists which yields maximum lift and' induced drag. 
This result is due to the fact that the effects of jet thickness and jet velocity are compen- 
sating. In other words, the thinner jet reduces the lift while the higher jet velocity 
increases the lift. With the jet thickness reduced to zero, power-off results would be 
produced independent of the jet-velocity ratio, as has been shown in reference 24 for the 
two-dimensional case. 

(b) When the jet Mach number is allowed to vary, as was done for the first and sec- 
ond curves, IV^ soon becomes too high as the thickness is reduced, that is, Mj g 1, for 
the linear subsonic theory to be valid. Therefore, the computation has not been carried 
out to find the thickness for the maximum lift and induced drag. 

(c) Accountir^ for Mach number differences resulted in increases in lift and induced 
drag with increases in Mj, which can be seen by comparing curves 1 and 3. This result 
can be explained by the reflection coefficient defined in reference 24: 
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where 3 = V 1 


Reference 24 showed that a more negative reflection coefficient will 


result in higher lift, and an increase in Mj tends to change Xj2 negatively. 

Effect of jet temperature .- The effect of jet temperature can be observed in figure 8 
by comparing curves 1 and 2 (solid and dashed lines) and noting that the density ratio 
( Pj j is equal to the ratio of the jet to the free -stream temperature, for the case where 
Pq = Pj. Increasing jet temperature tends to decrease the lift and induced drag, primarily 
because this increase would reduce the jet dynamic pressure. To see this effect, the 
ratio of free stream to jet dynamic pressures, T/i^, for small values of a can be 
reduced to 




1 + 1/1 + 


2Ct |T 


Ct|t / 

1 + — — + [i + 


2Ct -T 


P S 

1 2 
+ ■ > - + 


where equation (51) has been used. Equation (56) shows that Tp^ increases as T(= ^ 


increases. Increasing Tj would also affect Mj, but this would be a secondary effect 
compared to the dynamic pressure change. It is recalled that, for the present develop- 
ment, the ratio of specific heats in the jet is assumed equal to that in the outer flow' and 
thus is not a function of Tj, 

The lift reduction due to hot jet at a constant C-p has been measured experimen- 
tallv (ref. 8). 
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Effects of jet aspect ratio .- The jet aspect ratio is defined as the ratio of jet width 
to jet thickness. To investigate the effect of the jet center line is assumed fixed 
= 0.282^ and the jet span and the thickness are changed in such a way that the jet 
cross-sectional area is imchanged. The results at C-p = 0.2 are plotted in figure 9 for 
= 6.825 and 3.033. The jet-deflection angle was 10°. These results tend to con- 
firm the results of reference 10 which show that higher Alj will produce higher 
Cl r and Cp at a given a, where Cl p is the circulation lift coefficient. From 
the Cl l against Cjj ^ curve, however, the performance of the wing with either of 
these two jet nozzles is seen to be approximately the same. It should be remembered that 
higher ARj jet would yield higher scrubbing drag. 

Ef fects of jet def lection on cruise performance .- Again, the same cambered wing as 
used previously was employed with the jet exit either at 0.25c or at the leadir^ edge. As 
shown in figure 10, the jet-deflection effect is seen to be beneficial and the jet exit should 
be as close to the leading edge as possible. For the latter case, the jet will enhance the 
loading near the leading edge to produce higher leading-edge thrust. It will also produce a 
force component in the thrust direction due to the pressure loading at small a because 
of the positive camber slope there. At high a, this latter benefit disappears so that the 
jet tends to produce more induced drag than the power-off case. Since the jet-deflection 
effect is favorable on cruise performance, a thick wing which has high upper-surface 
slope at the trailing edge may be preferable to a thin wing from the aerodynamic point of 
view. 

Effects of jet interaction on the aerodynamic center .- At small values of the 

effect of power on the aerodynamic center is small, as has been found from the computa- 
tions at Crj, = 0.2. To determine the effect at high values of C,p, the configuration shown 
in figure 7 is again used, with the flap-chord ratio of 0.3, full-span flap deflection of 30°, 
and the jet exit at the leading edge. The aerodynamic center (a.c.) was computed by 
finding the ratio, -ACm/ACL, with ACm and ACl calculated between a = 1° and 5°. 
The jet-reaction effect is not included. As shown in figure 11, increasing jet thrust moves 
the a.c. forward, where the a.c. is expressed in terms of fraction of the root chord and 
is measured from the wing vertex. This conclusion is consistent with the experimental 
results (ref. 1) where the a.c. shift as high as 10 percent of the mean aerodynamic chord 
has been reported. Since the a.c. is affected by the jet- deflection effect, which varies with 
the flap angle, it follows that the a.c. shift must also be a function of flap angle. This 
effect is illustrated in figure 12 for C-p = 2.0, where Aa.c. represents the shift in a.c. 
location between power on and power off and is shown plotted with respect to 6j. The 
a.c. is shifted forward when the jet-deflection angle is increased. 

Effects of vertical shift of jet location .- The effect of jet vertical location on the 
aerodynamic characteristics is shown in figure 13, where the jet lower surface is taken 
to be on the upper surface of the wing and then is shifted upward one-half of a jet thickness. 
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It is asstimed that the jet -deflection effect is completely lost once the jet is shifted away 
from the wing surface. It is seen that at a given circulation lift the wing with the verti- 
cally shifted jet wovild have higher induced drag. This result is mainly due to the high 
rate of lift reduction associated with this vertical shift. 

To show how fast the aerodynamic characteristics would change as the jet is shifted 
upwards, the configuration of figure 5(a) without flap deflection and with the jet exit at the 
leadir^ edge was used. With Crp = 2.095 at a = 5°, the results of figure 14 show the 
rapid changes in AC^ and ACq ^ which are the net changes with respect to the wing- 
alone case. 

Figure 15 illustrates the effect jet vertical location has on span loadir^, where the 
left engine is inoperative. Shifting the jet vertically results in significant reductions in 

which leads to corresponding reductions in rolling- moment coefficient. Of course, 
the lift coefficient is also reduced. 

Effects of wing aspe ct r atio oi^lift cap ability .- The effects on lift a\:igmentation of 
changing wing aspect ratio from 8 to 4 are shown in figure 16 for a rectangular -wing plan- 
form. The aspect ratio is changed by varying the wing span with fixed chord length and 
jet geometry. The area of the wirg with ifl = 8 is used in computing Cl for all 
wii^s. The results were obtained at 0=5° and Sj = 30° = O.s) with the jet exit 

located at the leading edge. The results show that the reduction in circulation lift when 
the aspect ratio is reduced from 8 to 4 is 24 percent at C-p = 2.0. 

An analogous situation exists in reference 21 for a thin jet blowing at the trailing 
edge of a rectangular wing. For this case, the circulation lift for a wing of Al = 4.15 
with full-span blowing is about 31 percent less than that for a wing of Al = 8.3 with 
inboard half- span blowing. 


CONCLUDING REMARKS 

A theoretical method has been formulated which accounts for the wing -jet interaction 
and Mach number nonuniformity for the USB configurations. The theoretical results were 
in good agreement with experimental lift, induced-drag, and pitching-moment data for 
configurations with thin and thick jet exhausts. 

Additional theoretical results showed that (1) the lift and induced drag are reduced 
by increasing jet temperature and are increased by increasing jet Mach number; (2) reduc- 
ing jet aspect ratio with a constant jet cross-sectional area reduces the lift, induced drag, 
and pitching moment at a given ar^le of attack but with a minimal change in the induced 
drag polar curve; (3) the jet-deflection effect is beneficial at cruise and therefore a thick 
wing is preferable to a thin wing because of the greater upper -surface slope at the wing's 
trailing edge; (4) the aerodynamic center is shifted forward by power and/or jet-deflection 
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angle; (5) movii^ the jet away from the wii^ surface will change the lift and induced drag 
rapidly; and (6) reducing the wing span of a rectangular wing by half decreases the jet- 
circulation lift by only 24 percent at a jet thrust coefficient of 2.0. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
June 9, 1975 
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APPENDIX 


EXPRESSIONS FORINFLUENCE-COEFFiaENT MATRICES 


For the purpose of satisfying the boundary conditions, the continuous wing and jet 
vortex sheets are replaced by stepwise constant vortex distributions in the spanwise 
direction. The induced velocity vector due to a "bounded” vortex element of strength 
y dx' can be shown to be (ref. 17) 


where 



jg^ydx' axj fv 
|r xT|^\|b’| 



• i' 


a = (xi - x^T + ^yj - y^j + (z^ - z^k 
b = /x2-x)T + (y2-y)T + (z2-z)k 
a’ = ^Xj - x^i + ^^yj - y^j + jS^z^ - z^k 

b' = (x2 - x^i + - yp + j3('z2 - z^k > 

(X2 - Xj)T + (y^ - yj)T + (zg - z^)k 

= ("^2 ■ ^(^2 - yi )^ + ^(^2 - 

R = xi + yj + zk y 


(Al) 


(A2) 


The geometry is illustrated in the following sketch: 


For simplicity, define 


(xpypzi) 



(x2,y2,Z2) 



(A3) 
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Integrating equation (Al) over a vortex strip gives 


2 

^l(^) = It r(x‘)M/^x’,R) dx* 


(A4) 


where Xj, and x^^ are the end points of the vortex strip. Using the following 
transformation 


x’ = Xjg + i(x^ - Xj^)(l - cos 6) 


(A5) 


equation (A4) becomes 


5i(Ri) = — 5i^ sin e de 


^2,x _ X ) N 

k— ^ 5 Z 


k=l 


(A6) 


where the midpoint modified trapezoidal rule has been used to reduce the integral to a 
finite sum (ref. 17) and 


fl 2k - 1 

y a jr 

k 2N 


(k = 1, 2, . . N) (A7) 


Xi = 


2^ t 




{i= 1,2,. . .,N) (A8) 


In addition, the induced velocity due to the associated trailing vortices from end 
point 2 can be written as 



where 


R’ - R = (x* - x)T + (y' - y)T + (z’ - z)k 
r| = (x - x’)2 + I3\ - y’)2 + /B^Cz - z’)2 


(A9) 


(AlO) 

(All) 
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dFg = dXgT + dygT + dzgk 


(A12) 


Define 


G (x’,r)= y 
^(2) 


«> (r' - r) X dF, 


(A13) 


R 




For trailing vortices which are assumed to trail back on the wing plane, can be 

shown to be (ref. 15) 


^(2) ^ |32 


- k 


Zg - z 

A- 
1 1 

Xg-X 

(y - y^f * - ^ 2 )^]' 

(X - ^2f + - ^2f + 

V2 - y 

J. 

Xg - X 


_1 


1/2 


(y - Y^f + (z - zg'^ 


x 2 ^ 2 , x 2 „ 2 , 

X - Xgj + 13 fy - Y2) + i3 (z - Zg 


,ll/2 


(A14) 


Similarly, for trailing vortices emanating from end point 1, the induced velocity is 
given by (^noting that df j = -df2) 


qg \R)= ^ y(x’) dx’ G. (x’,R) 

47T 


'( 1 ) 


X, 


(A15) 


where 


pOO 

G. (x’,R = \ 
^(1) ^x’ 


^R’ - r) X dFj 


r; 


= -] 


Zi - Z 


( y - yi ) +( z-zi 


1- 


Xj - X 


(x - Xj)^ + ^^(y - y^)'“ + ^^(z - Zj) 


nV2 


+ k 


Yi - y 


Xj - X 


( y - yi ) 


X - x^f + / 3 ^(y - yj)^ + ( 3 ^(z - z^)' 


l/2j 

(A16) 
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G(x',R)= G (x’,R)+ G (x’,Rj 

( 1 ) ( 2 ) 


(A17) 


Then the induced velocity due to the trailing vortices from both end points is 


^2(^) = I? I ^ y(x’)G^(x’,R) dx' 


(A18) 


Again, applying the transformation of equation (A5) to equation (A18) 


(A19) 


is obtained. The total induced velocity, Yq, at a control point i due to a vortex strip 
is, therefore, the sum of q^^ and q 2 i 


n-fVii) y 

1/ 8N Ij 


k— 1 


(A20) 


It follows that the normal velocity induced at i due to a unit vortex strength at k, or the 
(i,k) element of the influence- coefficient matrix [^, can be written as 


(3 Ixf. - Xjj 




(A21) 


For the case where the unit normal vector n = k, it can be easily shown that 




(^1 - ^i)(y2 - yi) - (yi - yj)(^2 - ^i)fb' . t r .t 

? V |b'| |a"| , 


(A22) 
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Similarly, the (i,k) element of the matrix [s] can be written as 


ik 8N k 




(A24) 


For a = T, as being assumed in the present applications, it is known that 


G, ■ 1 = 0 


(A25) 


and 


M, • ej = 


_ (S'! - nX'a - "i) - (^1 - - >'i)fb' , ? ?.rA 


a X je' . 


b’ 


(A26) 


On the jet side panels, n would be j on the outboard panels and -j on the 
inboard panels. It follows that 


• n. = q: 


(*1 - =‘l)p2 - ^l) - (“l - ^i)(='2 - *l)/b' X' . ?■' 


a' xT'l^ 

'l 


b’ 


(A27) 


and 


• n. = ±. 


Z2 - z- 


yi - y2f + (^i - ^2? 


^ 1- 


X2 -x^ 


Zf-Zi 


- 1 


(Yi - Yif + (Zi - zi)' 


(Xi - X2f + /S2(y. - ygf + _ 22^ 

^1 - Xj 

(xi - Xj)^ + /3^(yi - Yif + - Zi)' 


lV2 


lV2 


(A28) 


where the upper sign is for the outboard side panel. Equations (A25) and (A26) are still 
valid for this case. 
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(a) Planform. 


ODA Experiment (ref. 4) 
Present theory 



a, deg 

(b) Lift curves. 

Figure 4.- Comparison of predicted lift curves with experiment 
for the USB configuration of reference 4. 6f = 0°, 5j = 12°, 
Mo = M, = 0. Linear dimensions are in meters. 



(a) Planform. 



(b) Lift curves. 

Figure 5. - Comparison of predicted lift curves with experiment for the 
USB configuration of reference 1, Mq = Mj = 0, 6f = 6j = 30°, 

JR = 7.79, ^ = 0.3, NACA 65J-412 airfoil. Linear dimensions are 
in meters. 
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Figure 8.- Effects of jet thickness, temperature, and Mach number nonuniformity 
on Cj__ and Cq ^ for planform of figure 7 with flat airfoil. Jet exit at 0.25c, 
Mq = 0.3, a = 5°, and C>p = 0.2. 
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Figure 11,- Power effect on aerodynamic- center location for configuration of 

Cf n 

figure 7. -^ = 0.3, 6f = 30”, and Mg = Mj = 0. (Jet-reaction effect not 
included. ) 



6., deg 

Figure 12.- Effect of jet deflection on aerodynamic- center location for 
configuration of figure 7. = 0.3, Cip = 2.0, and Mq = Mj = 0. 

(Jet- reaction effect not included.) 
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Figure 15.- Theoretical effects of jet vertical location on span loading and rolling 
moment for configuration of figure 5(a) with jet exit at leading edge and left 
engine inoperative, a = 5°, = 6f = 30®, and Ct = 2.095. 
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Figure 16.- Effects of wing aspect ratio on lift augmentation of rectangular planforms 
with USB jet. = 0.3, 6j = 6f = 30°, a = 5°, and Mq = Mj = 0. Jet exit at the 
leading edge, NACA 65^-412 airfoil. (Results are based on wing area of JR = 8.) 
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